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Abstract. The nonlinear memory effect is a slowly-growing, non-oscillatory 
contribution to the gravitational-wave amplitude. It originates from gravitational 
waves that are sourced by the previously emitted waves. In an ideal gravitational- 
wave interferometer a gravitational-wave with memory causes a permanent 
displacement of the test masses that persists after the wave has passed. 
Surprisingly, the nonlinear memory affects the signal amplitude starting at leading 
(Newtonian-quadrupole) order. Despite this fact, the nonlinear memory is not 
easily extracted from current numerical relativity simulations. After reviewing 
the linear and nonlinear memory I summarize some recent work, including: (1) 
computations of the memory contribution to the inspiral waveform amplitude 
(thus completing the waveform to third post-Newtonian order); (2) the first 
calculations of the nonlinear memory that include all phases of binary black 
hole coalescence (inspiral, merger, ringdown); and (3) realistic estimates of the 
detectability of the memory with LISA. 



1. Introduction 

We often think of gravitational-wave (GW) signals as having an oscillatory amplitude 
that starts small at early times, builds to some maximum, and then decays back to 
zero at late times. For example, this is the standard picture of a waveform from a 
coalescing compact-object binary. However, this picture is incomplete. In reality, all 
gravitational-wave sources possess some form of gravitational-wave memory. The GW 
signal from a 'source with memory' has the property that the late-time and early-time 
values of at least one of the GW polarizations differ from zero: 

A/i™T= lim /i+,x(<)- lim h+,^{t), (1) 

where t is time at the observer. 

When a GW without memory passes through a detector, it causes oscillatory 
deformations but eventually returns the detector to its initial state. After a GW with 
memory has passed through an idealized detector (one that is truly freely- falling), it 
causes a permanent deformation — leaving a 'memory' of the waves' passage. High- 
frequency detectors like bars or LIGO are rather insensitive to the memory from most 
sources because the detector response timescale is generally much shorter than the 
rise-time of typical memory signals (the characteristic time for the non-oscillatory 
piece of the GW signal to build up to its final value). A detector like LISA is better 
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able to detect the memory because of its good sensitivity in the low-frequency band 
where typical memory sources are stronger!^ Note also that bars and ground-based 
interferometers like LIGO are incapable of 'storing' a memory signal because they are 
not truly free: internal elastic forces push a bar back to its equilibrium shape, and 
magnets on the LIGO test masses (as well as its pendulum wires) push them back 
to their equilibrium positions. Because its proof-masses are truly freely-floating, a 
detector like LISA could maintain a permanent displacement. However, this late- 
time displacement caused by the GW memory is not directly observable without 
information on the prior state of the detector: the spacetime metric near a detector 
long after a GW with memory has passed is equivalent to flat spacetime in non- 
standard coordinates (see Sec. VD of [7] for details). Rather, it is the build-up of the 
memory (the difference in the metric between late and early times) that is observable. 

The memory effect has been known since the 1970's [51 IH HU] in its linear form. 
The linear memory generally arises in systems with unbound components: a binary 
on a hyperbolic orbit (two-body scattering) [TT], matter or neutrinos ejected from a 
supernova [121 [131 [14] , or gamma-ray burst jets [15]. In the 1990's a nonlinear form of 
memory was discovered independently by Blanchet & Damour [16] and Christodoulou 
[TT] . The nonlinear memory arises from the contribution of the emitted GWs to 
the changing quadrupole and higher mass moments. As discussed by Thorne [18j . the 
nonlinear memory can be described in terms of a linear memory in which the unbound 
masses are the individual radiated gravitons. This implies that nearly all GW sources 
are sources with memory (even if the component objects remain bound). 

Because the memory is a non-oscillatory effect with poor observational prospects 
for LIGO and other ground-based interferometers, it has been largely ignored by 
the GW community. But for the most important (or most studied) GW source — 
quasi-circular compact binaries — the nonlinear memory has quite a large contribution 
to the time-domain waveform amplitude: in a post-Newtonian (PN) expansion of 
the waveform polarizations, the memory effect enters at leading- (Newtonian)-order! 
That the memory enters at such low PN order is related to the fact that it is a 
hereditary effect — the memory amplitude at any retarded time depends on the entire 
past motion of the source (and not just on the source's instantaneous retarded-time 
configuration) . In addition, the nonlinear memory is a unique nonlinear effect because 
its non-oscillatory nature makes it distinctly visible in the waveform. For these reasons 
the memory should be studied further and its prospects for detection reassessed. 

In the rest of this article I will briefly review the linear and nonlinear memory 
effects. I will then summarize my previous work in three areas: (i) computing the post- 
Newtonian memory corrections to the GW polarizations; (ii) calculating the evolution 
and saturation of the memory during the merger and ringdown of coalescing black 
hole (BH) binaries; and (iii) estimating the ability of LISA to detect the memory from 
supermassive BH binary mergers. Except for the presentation of the linear memory 
from hyperbolic binaries in Sec. [3] this conference proceeding concisely summarizes 
the results of references [T] [19] as presented at the 8"^ Edoardo Amaldi Conference 
on Gravitational Waves. Readers are referred to those references for a more detailed 
exposition. 

§ It has recently been realized that pulsar timing arrays (T] [21 E] (which operate at much lower 
frequencies than LISA) could also be sensitive to the memory effect [IIISIE] , although the expected 
event rates from supermassive black hole binaries are small. 



The gravitational-wave memory effect 



3 



2. Preliminaries 



We first introduce some formalisni that will be useful for our discussion of the memory. 
The GW polarizations are conveniently decomposed in a sum over (Z,m) modes as 

CO / 

h+-ih^=J2J2 /i""-2r'"(e, (2) 

1=2 m.= -l 

where _2F'™ are spin- weighted spherical harmonics, and (0, $) indicate the direction 
from the source to the observer. In a multipolar expansion of the GW field, the modes 
/i'™ are related to the radiative mass ([/'"*) and current (1/'™) multipoles via 

/i'" - ^ [C/'™(T^) - iV"^{TR)] . (3) 

Here R is the distance from source to observer, Tr is retarded time, our units assume 
G = c = 1, and our overall sign and normalization depends on our polarization triad 
and other conventions (see [71 for details). The moments C/''" and y''" are constructed 
from their corresponding symmetric-trace-free (STF) tensors of rank-Z: 

U^"' ^ AiULyL^* and 1^''" = VlJ^I™ *, (4) 

where Ai and Bi are ^-dependent constants and y'j^ are the STF spherical harmonicsljj] 
The radiative moments C/'™ and F'™ that appear in the wavezone expansion of 
/i+.x are related to a family of source multipole moments {Iim, Jim, • ■ •)• These source 
moments are defined in terms of integrals over the stress-energy pseudotensor of the 
matter and gravitational fields of the source through a multipolar post-Minkowski 
iteration scheme (see [5^ for a review or Sec. II of [7] for a brief summary). For 
example, the radiative mass moments are related to the source mass moments via 



'V)dT + + 0(2.5PN).(5) 



Here ij^^ is the l^^ time derivative of the mass source moment , the integral term is a 
1.5PN order tail term, A4 is the mass monopole moment, ki is an Z-dependent constant. 
To is an arbitrary timescale that disappears in physical observables, jg 
the nonlinear memory term (discussed below), and 0(2.5PN) refers to several types 
of terms that enter at 2.5PN and higher orders. 



3. Linear memory 

As a simple example of the linear memory, let's consider the waveform from a 
hyperbolic binary. The leading-order multipolar contribution to the polarizations is 

h+-ih.^ ^ ^_2r'™(e,$). (6) 

m— — 2 ^ 

For a Keplerian binary in the x-y plane with relative orbital separation r{t), total 
mass M = rui + m2, reduced mass ratio ry = TOim2/M^, and orbital phase angle 
<f{t), the mass quadrupole is — {^(^''^/^V3)vMr{t)'^Y2jn{n/2,ip{t)). For Keplerian 

II These are related to the ordinary scalar spherical harmonics via y'™ = yj™n^, where ni^ = 
^21*^22 ■ ■ ■ ^ii ^ product of I unit radial vectors. On the quantities Ul, Vi, and y'j^, L represents 
I tensor indices (e.g., for I = 3, — > l^abc)- 
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orbits with semi-latus rectum p, eccentricity eq, and true anomaly v ~ ip ~ ujp {ujp 
sets the periastron direction on the x-axis), the orbital motion is described by 



P ... VpAi 

and V — (fi — — 5 — . (7j 



1 + eo cos V 

The resulting waveforms are given by ([6|) and 

T M~ 



Ao' = -8^/^77^eo(eo + cost;), (8) 



4±2 = -^/y'y^eT'"''*^ [1 - + (1 + eo cos«)(l + 2eoe±''')] (9) 

For < Co < 1 these waveforms are clearly oscillatory. But for a hyperbolic orbit (eo > 
1, with ojp — 0) the phase angle approaches ip- = — — arccos(— eg"^) at early times 
{t — > — oo), while at late times {t — >■ +oo) it approaches ~ v+ = arccos(— eg"^). 
This difference in the late and early time values of the orbital phase angle yields a 
corresponding difference in the derivatives of the mass multipoles, 

a41^ =0, (10a) 



AI^1>, = ±il6^y p eg ' 

resulting in a memory in the GW polarization amplitudes (see figure [Ij. Note that for 
a parabolic orbit (eo = 1) there is no memory since the orbital phase angle returns to 
its early-time value. As in the above example, the linear memory always arises from 
a change in the mass or current source moment derivatives, Alj^^ or AJ^^ for I > 2. 

We can also derive the linear memory for an unbound system by solving the 
linearized, harmonic gauge Einstein field equations (EFE) for the space-space piece of 
the metric perturbation hjk- Ohjk — —IQirTjk. Here Tjk is the stress-energy tensor 
of N gravitationally unbound particles with masses Ma and constant velocities va^ 
hjk is the trace-reversed metric perturbation, and □ is the flat-space wave operator. 
Solving this equation (via the Lienard-Wiechert solution) and projecting to trans verse- 
traceless (TT) gauge yields p^lTS]: 

N r . . tTT 



1-VA-N 



(11) 



Here TV points from the source to the observer and A means to take the difference 
between the late and early time values of the summation. In this formula the masses 
and velocities could refer to (i) the masses and velocities of the pieces of a disrupted 
binary, (ii) a gamma-ray-burst jet [15], or (iii) the individual radiated neutrinos 
[m [13] or pieces of ejected material in a supernova explosion (see, e.g., [HI dU [32] 
for numerical simulations of supernovae that show a memory effect.) 



4. Nonlinear (Christodoulou) memory 

The nonlinear memory [231 1161 [PT] (often referred to as the 'Christodoulou memory') 
arises from a contribution to the radiative mass multipole moments that is sourced 
by the energy-flux of the radiated GWs. Consider the relaxed EFE in harmonic gauge: 
Dh"^ = —16ttt"^, where t"^ depends on the matter stress-energy tensor T"^ , the 
Landau-Lifshitz pseudotensor i^L : ^^'^ other terms quadratic in h"^ [24j . Of the many 
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Figure 1. Examples of gravitational- wave signals with memory. The left plot 
shows the waveform modes 72m for a hyperbolic orbit with eccentricity eo = 2 
as a function of time [see ((Sjl and surrounding text]. Note the linear memory 
present in the imaginary part of 722 ■ The right plot shows the polarization 
for an equal-mass binary black hole coalescence with (blue/solid) and without 
(red/dashed) the nonlinear memory. The oscillatory piece of h+ was computed 
using an effective-one-body (EOB) model for the h22 mode. The memory piece 
was computed by substituting this mode into I I21I1 . See |19| for details. 



nonlinear terms in ; there is a piece that is proportional to the stress-energy tensor 
for GWs: Tj^ — -^^^^"-j^fe, where is the GW energy flux and rij is a unit 

radial vector. When applying the Green's function D"^ to the right-hand-side of the 
EFE, this piece yields the following correction term to the GW field [25] : 

J^t -|TT 

(12) 



5h]^ = ± I dt' 



R 



dt'dn' (1 - n' • TV) 

where is the retarded time. The time-integral in the above equation is what gives 
the memory its hereditary nature: the memory piece of the GW field for any value 
of Tr depends on the entire past history of the source. Thorne [18] has shown that 
the nonlinear memory (jl2l) can be described in terms of the linear memory (jlip if the 
unbound objects in the system are taken to be the individual radiated gravitons with 
energies Ea — M^/(l — w^)"^^^ and velocities v\ = cn\. 

When the GW field is decomposed into modes as in ([2]), the nonlinear memory 
can be shown to yield a correction to the radiative mass multipole moments ([5]) that 
enters at 2.5PN and higher orders [l6l [7]: 



^(nonlinmom) _ / 2)! f f dE'gw 



The radiative current moments V/„j do not have a nonlinear memory contribution. 
Note that the energy flux is itself defined in terms of the waveform modes via 

^^ = |l(Ai + ii>^|l S <^--'"- )-.^'''"' (")-.>■'"■""•(«). (") 

/' .m' ,m" 

This means that the memory is calculated in an iterative fashion: those him modes 
that contain the memory are negligible in the computation of the energy flux. 

The memory (linear and nonlinear) should not be mistaken as a change in the 
monopolar-piece of the 1/R expansion of the metric. Rather it is a change in the 
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quadrupolar (and higher-order) pieces of the l/i?-spatial-part of the TT-projection of 
the metric. It is a purely GW effect, and is not directly connected with the change in 
the 'Coulomb part' of the metric or the 'mass loss of the source' (spherically symmetric 
mass loss produces no GWs!) except indirectly through the changing mass's effect on 
the quadrupole and higher-order multipole moments. 



5. Memory contribution to the post-Newtonian waveform of 
quasi-circular, inspiralling binaries 

To compute the nonlinear memory contribution to the waveform for quasi-circular 
binaries, one begins by substituting (fT4|) into ([T3| and evaluating the angular integral 
in (jl3p . Next one substitutes explicit expressions for the him modes into the result. 
These modes are proportional to him oc e"""'''^*'' and are given explicitly (to 3PN 
order and for to 7^ 0) in [26]. The resulting integrand of the time- integral in ([T^ 
will be either oscillatory or non-oscillatory. The non- vanishing, non-oscillatory modes 
turn out to be the ^ifith Z-even. (The to 7^ modes are oscillatory — they 

do not contribute to the memory, but they affect the waveform at 2.5PN and higher 
orders.) These consist of time-integrals of polynomials in a; = (Afw)^/'^ 

(where w is the orbital angular frequency), with each term of the form 

[x{t)X'dt^ / —dx. (15) 



After performing the change of variable indicated in (jlSp , evaluation of the resulting 
integrals requires a model for the frequency evolution of the binary. The adiabatic 
evolution of the frequency (or x) is easily derived from energy balance (GW luminosity 
equals the rate of change of the orbital energy, £gw = —E) and the relation 
X = — £Gw/(d£'/da;). The resulting frequency evolution (known to 3.5PN order), 

+ 0{x) + 0(x3/2) + . . . + 0(2:^/2) ^ o(4PN)l , (16) 



dt 5 



can be used to evaluate integrals of the form ([T5|). For details see [?]■ 

When the above procedure is carried out, the resulting memory contributions 
only affect the -I- waveform polarization^ and have the form 

^(,,em)^ 277Mxg^„/2^|„/2,,,e„.)^ whcrC (17) 

Tl = 

^''"''"^^ ^4(17 + 4), (18«) 

^(0.5,mcm) ^ q ^ _^(1.5,mcm)^ ^^3^^ 



^(l.mcm) 



354 241 62 059 , 4195 



e " coo 1 00*^0 



2 064 384 1032192 " 688128 
15 607 9373 ^ ^15_ 4 
73 728 ^ 36 864^® 8192^® ' ^ 



(18c) 



^ This is a consequence of our choice for the polarization triad. A rotation of the triad would cause 
memory in both polarizations. 
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(2,1 



3 968 456 539 570 408173 



9 364 045 824 
7 169 749 



4 682 022 912 
13 220 477 



Co 



18 579 456 18 579 456 



122166 887 
^ 3121348 608 
1345 405 
6193152*^ 



75 601 
15 925 248 



4 



25115 
884 736' 



10 097 



5179 



147456 36 864 



44 765 
147456 



4 + 



3395 
73 728 



H 



(2.5,1 



57r 



21504 



(1 - 477)s|, (509 + 472c| + 394) 



H 



(3, mem) 



Equation (4.6g) of [7], 



(18rf) 
(18e) 

m 



where cq = cosO and sq = sinO. 

The OPN memory contribution ()18ap was first computed in [25]; the 0.5PN term 
was computed in 26 . All the other PN corrections listed above are new and derived 
in [71 . Combined with the oscillatory waveform pieces computed in [26] , these results 
complete the waveform polarizations at 3PN order. 

Note in particular that at Newtonian (OPN) order, the total waveform is given by 



1 



-(1 + 4) cos2((^ + $) + ^4(17 + c|) + 0(xi/2) 



^(OPN) ^ ^rjM^ 
R 

This shows that even though the memory originates in a high-PN-order effect, it 
enters the waveform at the same order as the standard, leading-order 'quadrupole' 
piece. This fact is not widely appreciated. 



-2ce sin2((^-$) + 0(a;i/2) 



(19a) 



(19&) 



6. Computing the memory in binary black hole mergers 

6.1. Limitations of numerical relativity in modeling the memory 

The PN waveforms computed above allow us to model the slow growth of the 
memory during the inspiral phase of coalescence. To determine the rest of the 
memory's evolution and its eventual saturation value, we need to model the merger 
and ringdown phases. Numerical relativity (NR) provides the only accurate way to do 
this. Unfortunately, the memory is not a quantity that is easily computed by current 
NR simulations. Indeed, one may have noticed that none of the waveforms produced 
by recent NR simulations show any sign of memory. 

There are several reasons for this. First, for quasi-circular inspiralling binaries 
the memory is not present in the commonly plotted I — m = 2 mode. The memory 
is only present in the m — Q modes of the waveform. More importantly, while the 
memory in the metric perturbation and its /i;o modes can be comparable in size to the 
leading-order oscillatory /i2m modes [see e.g., (|19a|) ]. the memory contribution to 'I'4 
(the Weyl scalar that NR simulations usually use to extract GWs) is suppressed by 
several orders-of- magnitude relative to the leading-order oscillatory piece of ^1*4. For 
example at large i?, ^'4 — )■ — i/ix , and a decomposition of ^'4 similar to ([2]) implies 
that its modes tpim are related to the metric perturbation modes via tpim = him- 
During the inspiral the metric modes be can decomposed into an amplitude Aim{t) 
that slowly evolves on a radiation-reaction timescale T„ oc {M / vf){Mu})~^/^ , and a 
phase that oscillates on an orbital timescale Torb oc 1/w: him — Aim{t)e'~^™'^*- . Then 
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for oscillatory (m 7^ 0) modes il^im is snialler by two factors of the orbital time, 
t/j/m ^ uP'him, but for memory (m = 0) modes is smaller by two factors of the 
radiation- reaction time, ■0/0 ~ if' /M'^i^MuY^/'^hiQ. This means that the tfif^ memory 
modes are smaller than the oscillatory 'il'im modes by 5PN orders. The memory modes 
of \l/4 are therefore a small correction superimposed on the much larger oscillatory 
pieces of ^'4 and are thus difficult to resolve in a numerical simulation. 

Aside from this problem, there is also the issue of computing the two integration 
constants needed to go from to him- Berti et. al [27] showed that choosing these 
constants incorrectly leads to an 'artificial memory' that shows up in all {I, m) modes. 
One approach to computing these constants correctly is to match the NR waveform 
modes onto the PN waveform modes. 

Some of the above problems can be ameliorated by using a wave-extraction 
formalism that directly computes the metric waveform without the need to compute ^'4 
and integrate twice. However another problem arises that also affects metric-extraction 
schemes: the sensitivity of the nonlinear memory to the distant past requires large 
initial binary separations (> 50Af) to compute the hiQ modes accurately. See Sec. VC 
of [3 for further discussion of these issues. 

6.2. Analytic modeling of the full coalescence 

Since NR simulations of the coalescence memory are not yet available, we need to 
investigate other approaches. A first attempt by Kennefick [28] took the leading- 
order inspiral memory [ (jlOap . which grows as /j™°™''"''p q<; _ where tc is the 
coalescence time] and assumed that it abruptly stops growing at a time corresponding 
to a specified (and somewhat arbitrary) orbital separation (roughly the last-stable- 
orbit). This approach truncates the memory too early and does not incorporate its 
smooth build-up during the merger and ringdown. 

To improve upon Kennefick's analysis we begin by introducing the so-called 
minimal waveform model (MWM). The purpose of the MWM is to use information 
from NR and the effective-one-body (EOB) [29^ approach to construct a very simple, 
fully-analytic waveform that qualitatively contains the basic features of a binary BH 
coalescence waveform. The MWM simply consists of matching the leading-order 
inspiral waveform to a ringdown waveform. For example, the g"^ time-derivative of 
the (2, 2) mode is approximated as 



where u: = ^ = (M/r^y/^, r = r^il - (t - t™)/T„]i/4, r„ = (5/256)(M/r,)(r„/M)4, 



o'lmn = iwimn + are the BH quasi-normal mode frequencies and damping times 
[50] (these depend on the final BH mass and spin which are determined from NR 
simulations), and tm is a 'matching time' at which r = r^. The matching radius r™ is 
an adjustable parameter that determines the peak of the waveform amplitude; I choose 
the value r,„ — 3M corresponding to the Schwarzschild light-ring. The coefficients 
Aimn are determined by matching (pil)) at t = for < q < rimax- 

For quasi-circular orbits, the leading-order piece of the memory is sourced by 
the /122 oscillatory waveform mode. To compute its contribution to the memory. 




for t < t 



for t > t, 



(20) 
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substitute the / = 2 expansion of the energy flux (HH) into ^ to yield 

and the corresponding hiQ modes. The resuhing memory contribution to is 

hf-^^^sUn + cl) fjh22\'dt. (21) 

This reduces to a simple analytic expression when the MWM (1^0]) is substituted 
(see [12] for details). The effect of including higher-order waveform modes in ((Ti)) 
is currently under investigation. Because the memory (j2ip depends only on the 
amplitude of the /122 mode and not its phase, the MWM provides a reasonably accurate 
description of the memory. As an alternative to this simple analytic model, one can 
use an EOB description of the /122 mode and substitute this into (PT|) . This approach 
is more accurate because the EOB waveforms are calibrated to NR simulations. Both 
methods were implemented and compared in |l9j . Figure [1] shows the memory resulting 
from the EOB-based model. 



7. Detecting the nonlinear memory 

Although the memory has most of its power at near-zero frequency, the build-up of 
the memory distributes power in a range of frequencies. As a simple estimate of the 
memory's detectability with GW interferometers, we compute the sky-averaged rms 
signal-to- noise ratio (SNR) for a detector with noise spectral density 5„(/), 

ni/2 



SNR: 



hlif) df 



hlif) f 



(22) 



where the sky-averaged rms noise amplitude is hn{f) = afSn{f) {a = 5 for LIGO- 
like interferometers, 20/3 for LISA), and the characteristic memory amplitude is 

Kif) = 2(1 + z)f{\ur'^\ii + z)/]ni/^u^z,,/(i+.), (23) 

where is the Fourier transform of Dl{z) is the luminosity distance, 

and the angle brackets mean to average over sky-position and polarization angles. For 
the MWM the Fourier transform for / > can be expressed analytically as [T9] 

r(mcm,MWM)/ „s V^I 2 r. , 2 \ '^^ J^'^^'^Fl O -f TTf^ -7 I A O M 

K 'if) = 3847ri^ ^e(17 + ce)^< — — [1 - 27ri/T„.C/(l, 7/4, 27ri/T„.)J 



1 '''22ri'^22n'^22n^22n' 



nTto '^^'^ - (^22„ + ^^2 



(24) 



where U is Kummer's confluent hypergeometric function of the second kind and J- is 
a 'fudge factor' that matches the MWM and EOB memory models at late times. 

Evaluating the SNR, one finds that initial LIGO can only detect the memory 
from stellar-mass BH mergers within the Local Group. Advanced LIGO will see the 
memory to about 10 times that distance (< 20 Mpc). LISA, however, will see the 
memory from supermassive BH mergers out to rcdshifts z ^ 2; see Figure [H 

8. Conclusions 



Gravitational waves with memory are potential sources for detectors with good low- 
frequency sensitivity (LISA and pulsar-timing arrays) and are worth further study. 
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Figure 2. Dctcctability of the memory with LISA. The left plot shows the LISA 
noise amplitude hn (soUd/black) and the characteristic amplitudes he for the 
inspiral waves (solid/orange) and the nonlinear memory for a black hole binary 
with the indicated parameters. The short-dashed/blue curve uses the MWM 
II24I I to compute the memory; the long-dashed/red curve uses a truncated-inspiral 
model based on 1281 (see |19l for details). The right plot shows the angle- averaged 
SNR of the nonlinear memory signal for equal-mass LISA binaries as a function 
of the total binary (source-frame) mass M and redshift z. 



The nonlinear memory is particularly interesting because (i) it arises from the 
gravitational "waves produced by gravitational waves, (ii) it affects the "waveform at 
leading-order, and (iii) it imparts a unique and visually apparent signature to the 
"waveform. Although our kno"wledge of the oscillatory pieces of the waveform has 
greatly increased in recent years (thanks to the work of many in the PN and NR 
communities), our understanding of the memory has not. Here I summarized recent 
attempts to fill in the gaps in our understanding. This included determining the 
missing pieces of the 3PN quasi-circular inspiral waveform that arise from the memory. 
It also included a first study of the build up and saturation value of the memory during 
the merger and ringdown phases. As discussed above, numerical relativity simulations 
cannot yet easily compute the nonlinear memory. This study therefore serves as 
an example of how post-Newtonian and effective-one-body methods, combined with 
crucial input from numerical relativity simulations, allows us to compute quantities 
that are not directly computable with analytical or numerical techniques alone. The 
prospects for detecting the nonlinear memory were also investigated and found to be 
somewhat poor for advanced LIGO but good for LISA. Motivated by this work, several 
other studies [H O [B] have recently investigated the detectability of the memory with 
pulsar timing arrays. Work is in progress to improve upon the initial estimates of the 
merger/ringdown memory discussed here, extend some results to eccentric orbits, and 
perform a more detailed study of the memory's detectability with LISA and future 
ground-based detectors. 
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